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We carryout a comparative study of spin distributions defined over the sphere for bipartite 
quantum spin assemblies. We analyse Einstein-Podolsky-Rosen-Bohm (EPRB) spin correlations in 
a spin-s singlet state using these distributions. We observe that in the classical limit of s — > oo, 
EPRB spin distributions turn out to be delta functions, thus reflecting the perfect anticorrelation 
property of two spin vectors associated with a spin-s singlet state. 

^ ■ I. INTRODUCTION 

o: 

^^ ' Distribution function description of quantum mechanics owes its origin to the work of Wigner |1| . This description 

offers a framework in which quantum phenomena can he discussed using as much classical language as allowed. It 
O i' appeals naturally to one 's intuition and can often provide useful physical insight that can not be easily gained through 
1^ I other approaches 0. Irrespective of a few shortcomings, Hke for eg., appearance of negative probabihties, advantage 
of using the distribution function approach is in the fact that it involves classical functions as opposed to operators. 
fSJ , Besides the phase space distribution functions, where position q and momentum p are the statistical variables, there 
CN ■ have been several attempts to construct the spin distribution functions([3]-[6]). Infact, Wigner like distribution 
functions have been used ([5], |7|) to discuss EPRB spin correlations in a spin-s singlet state. In this formulation, 
" ' , the spin correlations are cast in a structurally similar form to those of local hidden variable models but with negative 
J^ ■ distribution functions. 

Qv , In this paper we extend the spin distribution functions defined over the sphere viz., the P— , Q—, and F— functions 

^^ ' ([4] I [6]) for bipartite spin systems characterized by a mixed state spin density matrix 5. In Section II we discuss 
^\ ' the expansion of the spin density matrix in terms of Fano statistical tensor parameters ||8| . In Section III we express 
^^ . the P— , Q—, and F— spin distributions in terms of Fano statistical tensor parameters. In Section IV we extend 
^^ ' these distributions to bipartite spin systems and study the EPRB spin correlations using them. We also show that 
^^ , distribution functions characterizing EPRB spin-s singlet state approach delta function in the classical limit of s — > cx). 
^ ' 

^ . II. FANO STATISTICAL TENSOR PARAMETERS 

^ ' A set of (2s -t- 1)^ spherical tensor operators {'r^(S), k — 0, 1, 2, . . . , 2s and q = — fc, — fc + 1, . . . , k} constructed 
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out of the spin operator S through 









2n4^(2s-fc)!(2s + l)- ^2.1) 



{2s + k + l) 



provides a linearly independent, orthonormal basis and any operator, acting on the (2s + 1) dimensional Hilbert space 
of a spin-s assembly, could be resolved into its irreducible components in terms of this basis. Here, Ykq{9, (j)) denote 
spherical harmonic functions and the normalisation constants Afsk are chosen so as to be consistent with the Madison 
convention fR!ol: 



(sm'|fg^(S) \sm) = V2fc + 1 c{sks;mqm'), (2.2) 

where c{sks; mqm') denote Clebsch-Gordan coefficients. From the property 

s 

y^ c{sks;mqm) ^ {2s + l)SkQSqQ (2.3) 
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of the Clebsch-Gordan coefficients, it could be easily verified that 'r^(S) are traceless for all non zero values of k and 
Tq (S) is a (2s + 1) X (2s + f ) dimensional unit matrix. 

By construction, t^(S) arc irreducible under rotations, with the transformation property: 

rotation r ^ n rotated 
f,^(S) _^ [t^{S)\ =i?(a,/3,7)r,^(S)i?t(a,/?,7) 

k 

- E D'^,',{c.,(3,j)f^,$), (2.4) 

q' — — k 

where D'' denotes (2k + 1) dimensional irreducible representation of rotations O] and a, f3, 7 are the Euler angles 
of rotation. Hermiticity property of f^(S) is expressed through, 

f,^g)t = (-irf'!^(g). (2.5) 

The orthonormality property, 

Tr (t^,$) f,V(g)t) = (2s + 1) 4fc'^„', (2.6) 

could be readily realised by making use of Eq.(2.2), (2.5) and the properties, 

c(siS2s; 17111712171) = (-l)"!-'"! J j^;^ c(siss2; mi - m - 7712), 
c(siS2s; mii7i2m) = (— i)«i+s2-s c(siS2s; — rrii — m2 — m), 

c{siS2S;i7li17l2l7l) — (— l)'*^^*^^* c(s2SiS; TO2TO1TO) , (2-7) 

together with the orthogonality 

Sl S2 

y^ y^ c(siS2S;mim2m)c(siS2s';mim2m') = (5^5' (5™^', (2.8) 

mi = — Sl 7n2 = — S2 

of the Clebsch-Gordan coefhcients. 

Any arbitrary operator A{S) can be expressed in terms of the irreducible tensor operators t^(S) as 

^ = EE ^'(S)t4 (2.9) 

k=0 q=-k 

where the spherical componets a^ of the operator O are given by 

a^ = Tr(6f,^-(S)). (2.10) 

Specifically, the spin density matrix p, which characterises a spin-s assembly, has the resolution: 

^ = ^EE^'(s>t^-, (2.11) 

k=Oq=~-k 

and i^' — (f^(S)) — Tr (pf^(S)) are the well-known Fano statistical tensor parameters [8]. 

Depending on the requirements of the spin density matrix p, the Fano statistical tensor parameters t^ satisfy the 
following properties: 

• Normalization: Tr (p) = 1 =^ t" = 1, 

• Hermiticity: p^ = p =^ t^* = (— 1)"?^^^, 

• Property under rotation: 



rotated 

k 



k, q 



k, q 

q' — — k 

For an entangled spin assembly containing a pair of subsystems with spins si and S2, the above discussion can be 
readily generalised( |12[ , |l3]), and the spin density matrix characterising such a composite bipartite system is given 
by, 

^- = (s^tttiWtt) £ i: £ i: 4t^*(^^^(so«^^(s;)), (2.12) 

^ '^ ^ ki=0qi = -kik2=0q2 = -k2 

where tg^^^ denote coupled Fano statistical tensor parameters, Si , S2 denote the spin operators of the subsystems 1 
and 2 respectively. The reduced or the subsystem density matrices are obtained by taking the partial traces: 

Pi = Tr^ (/5i2) = ^^-^ f: E rj^^H^iYt'S (2.13) 

and 

^^ ^ ^""^ ^^^^) ^ 72^ £ £ -,'MS>C- (2-14) 

*> 2 ) k2=0q2=-k2 

where we have made use of Eq.(2.2) and (2.3). The system is said to be entangled iff 

P12 ^p,^p,=^ e^ ^ 4^0 Q . (2.15) 



III. P- Q- AND F DISTRIBUTIONS IN TERMS OF FANO STATISTICAL TENSOR PARAMETERS. 

The coherent state description of electromagnetic fields has proved to be successful in providing an insight into the 
relationship between semiclassical and quantum theories of light ||14| .The diagonal coherent state representation or 
P-representation for the density matrix, 

p= [ da P (a) \a){a\, (3.1) 

has proven to be very useful in bringing many of the results of quantum electrodynamics into forms similar to those 
of classical theory and the expectation values of any quantum operators A could be realised as classical averages in 
terms of the weight function P{a) through 

i\=Tr pA = f daP{a)A{a), (3.2) 

where A{a) is the classical function corresponding to the operator A. Arecchi et.aL,[4] introduced the analogue of the 
diagonal coherent state representation for the spin density matrix p via the relation 

p = ldnp{e,c^)\e<f>){ecp\, (3.3) 

where dfl ~ sm9d6d(p and and 16*0) represents the spin coherent state(SCS) or Bloch state, defined as a rotated 
maximum 'down' spin state |s — s) : 



,(rS+-r*S_) ls_.s\ = 



E 



g,, ^^ . ,.-j |_^ _ ^^^ ^ _S(0- 7r,6',7r - (/)) \s - s) 

s-\-'rn 



2s 

S + ?7l 



COS- I I sm- I 6-*'"+™)'^ Ism) 



(3.4) 



Here, t = ^6 e "^ and S± are the spin ladder operators; R{(j) — tt, 0, tt — (/>) denotes rotation through Euler angles 
(j) — TT, 6 and tt — (j). From the normalisation condition Tr [p] = 1 we have 



dnp{e,i 



(3.5) 



i.e., the weight function P{9, (j)) in the diagonal spin coherent state representation is a normalised function. With the 
help of P{9, 4>), the quantum expectation value of any arbitrary spin observable A is given by the classical average 



A 



dnp{d,4>)A{0,(i>), 



where 



A{d,4>)^(0(t>\A 



(3.6) 



(3.7) 



is a classical function corresponding to the quantum mechanical operator A. For example corresponding to the spin 
operator S = (Si, S2, S3) we obtain 



S(0,0) = (00|S 



s n(7r — 9, 1 



(3.8) 



where n(7r — 9,(j)) = (sinScosc/), sin 6* sin ^, —cos 6*) is a unit vector defined on the Bloch sphere. Thus in the P- 
representation of spin, one can visualise the quantum expectation value of spin as a classical statistical average of 
random orientations of angular directions n(7r — 0, 0). 

Observe that the Fano statistical tensors t^ can be expressed as the classical averages 



i, = Tr 



pf^is) - dnp{0,cf,) {0cb\f^{s)m. 



(3.9) 



The expectation values of the irreducible tensor operators t^(S) in the spin coherent states \0(j)) can be simplified as 



follows: 



ffc (S) \0^) ^ (s - s| i?t (0 _ TT, 61, TT - </>) f^ (S) i?(</. - tt, e, tt - 0) |s - s) 



J2 {s-s\f^,{S)\.s-.s) D^^;, 

q' — — k 



TT, W, TT 



= ^ c{sks;-sO- s)V2kTlSq> oDg*,{(j)~'Tr,0,n 

q'^ — k 

{-!)'' c{sks;s0s)V^Ykq{9, (I) ~tt), 



(3.10) 



where we have made use of the transformation property of t^(S) under rotations. Utilising the explicit expression for 
the Clebsch-Gordan coefficient c{sks; sOs) Q given by. 



c(s/cs;sOs) = (2s)! 



(2s + 1) 



we obtain. 



{S)\0^)=V4^{-l)''+'^{2s)\ 



{2s - ky.{2s + k + ly. 
(2s + 1) 



(2s-/fc)!(2s + fc+l)! 



Ykg{0,^), 



(3.11) 



(3.12) 



where we have used the symmetry[ll] Ykq{0, 0— tt) ~ (^1)'' Ykq{0, 0), of the spherical harmonic functions. Substituting 
Eq.(3.12) in Eq.(3.9), we observe that P{0, </>) must be of the form 



2s k 
I 

a/zItT 



^E E {~lf^'Vsut\Y*^{e.<i^). (3.13) 

V47r ,^ — , 

in order to reproduce the Fano statistical tensors t^ in this representation. The coefficients Vsk are given by, 



_ 1 (2.-fc)!(2. + fc + l)! 

^^^-(2^V I2ITT) • ^^-^'^ 

Another useful distribution function which can be derived using SCS is the positive, normalised Q-function[2, 5]: 

Q(^'<^) = ^^^WI/5|^<^>- (3.15) 



Using Eq.(2.11) and (3.10) we obtain, 

2s k 
An \"firi"r/ ^^ 



Qi(^^^) = ^^^ mpm - i^E E *r mr^is)m 



k=Oq=-k 

it E ^^rMf->^^.(M)(2.)! Jj^^_,^,(2^^,^^), 

K/ — U Q — " — fZ I 



(2s + 1) 



2s k 

47r 

k=Oq=-k 

where 



i=E E (-1)'" 2" »5^>,(«'*) P-w) 



It could be readily verified, using the orthonorniality property of the spherical harmonics [11], 

''dn Ykg{e, 0) Yk'q' (9, (j,) = 6k, k' 5q, q. , (3.18) 



that the P— and Q— functions [|6[ given by Eq.(3.13) and Eq.(3.16) respectively, are normalised: 

jdnp{e,4>) = -^E E (-1)'^' ^^'^- tqJdnYC^ie,(t^) 

k—0 q— — k 
. 2s k 

^T.T. (-1)'^' ^^'^ 4 ^^ ^i^' oSq,o = l (3.19) 

^''^ k=Oq=-k 



and 



2s k 



dn Q{d, cj,)^ 4^ E E (-1)'^' Q^^ < Jd^ ^4(^' ^) 

^ k—Q q—~k 

-. 2s k 

= ^^T.T. (-1)'^' 2^^ ^9 ^^ '^fc. oSq,o = l (3.20) 

Also, the reality of the sums \^ (~1)^ tq ^fc*g(^i<^) ensures that P{6,(j)) and Q{6,(j)) are real. 

One can also construct a spin distribution function over the sphere using the characteristic function approach[6]. 
In this method, a classical distribution function F{X) with X = (ATi, X2, X^) as the associated random variables, is 
constructed by taking the Fourier inverse of its characteristic function (/>(/), which is the expectation value of e* 
i.e., 



(f>{I) = E[e''-^) = d''XF{X)e''-^. (3.21) 



Here, E {. . .) denotes the expectation value. It can be readily seen that the distribution function F{X) is obtained as 
Fourier inverse of the characteristic function (j){I): 

Fix) = j^JJJ d?Ie-^^-^m. (3.22) 

Making use of the well-known expansion 

oo k 

e-^"*-^ = 47r^ ^ z'^ j^{IX) ^^(0,0) Ffc,(07,07), (3.23) 

where, jk denotes the spherical Bessel function and the spherical polar co-ordinates of / and X are denoted, re- 
spectively, by (/, 6*7, (/)/), (X, 0, 0). One can express the characteristic function 0(/) in terms of the spherical 
moments 

,i\=E[ju{IX)Y^,{e,cj,)\, (3.24) 

as, 

oo k 

m = 47r^ ^ {i)^ YC^iOi, <i>i) fi"^^. (3.25) 

k=Oq=-k 

In the case of quantum spin assemblies, X corresponds to spin operators S = (5*1, 5*2, 5*3) and the classical average 
E {. . .) is replaced by the quantum mechanical average Tr (p. . .). Moreover, a spin distribution is expected to auto- 
matically reflect the constancy of the squared angular momentum i.e., S'^ = S^ + S2 + S^ = s{s + 1). On imposing 
the condition X^ + X^ + X^ = R"^ = s{s + 1) and introducing solid harmonics ykq{X) = R^Ykq{9, 4>) in Eq.(3.24) we 
obtain, 

M^- = jkilR) R-'' E (ykqiX)) , (3.26) 



where jk{IR) R ^ has been taken outside the expectation value since R = \/s{s + 1) is a constant. We can now use 



the correspondence rule |17 



y^.i^)^l^^m (3.27) 



so that 

for k <2s and E {ykq{X)\ — for k >2s, in accordance with the Wigner-Eckart theorem[ll]. Thus the characteristic 
function of Eq.(3.25) takes the form, 

2s k 

m = 47r^ ^ (,)fc _^i?-fe^,(/i?)r^;(0,,0,)ife. (3.29) 

k=Oq=-k ■ '^'' 

The Fourier transform of 4>{I) can now be readily obtained by making use of the rcsult[15], 

j J J I^dl sm0i d9i d^i e-^'-'^jkilR) n;(^/, h) = -j^i-^f S{R - X) Y,*^{e, cj,) (3.30) 



and we obtain the spin distribution function over the surface of the sphere of radius R = \/s{s + 1) 



Fix) 



OO fTT f27r 



(2^)3 ^ J, J, 

2s k 



rdl sm9ideid(j)i e-"^'-^ (l){I) 



5iR-x)j2Y: fcrirr^"'"*'^'^*^^^''^)- 



(3.31) 



The normalised distribution F( 



F( 



k=Oq=-k 

of angular co-ordinates 0, (j) is obtained through the relation, 

2s k 

- kq \ 






fe=0 q=-k 



where 



•^"'^ ~ 2k 



(2s + fc + l)! 



47r(2s - fc)!(2s + l){s(s + 1)}*^ 



(3.32) 



(3.33) 



It could be easily verified that F{9, (p) is a real, normalised function of 9, (j). All the three functions P{9, (f)) 

and F{9, (f>) have the properties of statistical distribution functions in the sense that, they are normalised to unity 

and they yield correct expectation values for quantum mechanical spin observables, based on the correspondence 

rules associating t^(S) to spherical harmonic functions Ykq{9,(f)), each with different weight factors, given explicitly 
through, 

• P-representation: 



r'{S) 



4. (,^^^^), (^l)-«n,(^,. 



(2s + 1) 



^(2s-fc)!(2s + fc+l)!' 



(3.34) 



• Q-representation: 



r^iS) 



/l^W^^^'^"'^'^'^ + ' + '^''^^'^''^^' 



(3.35) 



• F-representation: 



r'iS) 



47r(2s + l)(2s-fc)!{s(s + l)}* 
(2s + fc+l)! 



Y, 



kq[ 



(3.36) 



IV. P-, Q-, AND F- FUNCTIONS FOR BIPARTITE QUANTUM SYSTEMS 

The P and Q distribution functions can be readily extendend for bipartite spin assemblies as follows: 

p = I Jdni dn2 P{9i,^i;92, 02) l^i^i; 02, 02) (01, 0i; 02, 02| (4.1) 

and 

g(0i,0i;02,02) = ^^^i^^^|^^^(0i,0i;02,02|/5|0i,0i;02,02) (4.2) 

where |0i, 0i; ^2, ^2) are product spin coherent states with (^i, (jji), (^2, 02 ) representing angular directions on respec- 
tive Bloch spheres corresponding to systems with spin si and 82- The results of section III could be easily generalised 
to the case of a pair of systems and we obtain 

2si fei 2s2 k2 

p(0i,0i;02,02)--E E E E (-1)'^^'^ (-1)'^^'^ ^^^^^ ^^^^^ 



fci=Oqi= — fci k2=0 q2 = — k2 



>^t'^Y,\,S(^u<f>i)Y*^^^{02,h) (4.3) 

and 

2si ki 2s2 k^ 

ki=0 (ji = — fei k2=0 q2 = — k2 

^t'l^^yk*,JS,,MYC^je,,cl^,) (4.4) 

Let us now consider a EPRB spin-s singlet. The spin density operator in this case is p = |(ss)00) ((ss)00|, the 
matrix elements of which are given by 

{sm[; sm'2\p\smi; sm2) = —-- — — -— S^' -m' ^mi -ma, (4-5) 

(ZS + 1) 12 

in the (2s + 1) x (2s + 1) dimensional direct product spin space of particles 1 and 2. The coupled Fano statistical 
parameters t'^lq^ characterising EPRB spin correlations are given by, 

4^',^= Tr[pf,'^/(s"i)xf^^^(S2)] 

E (2s + 1) ^"^[ -m'2 ^-mi ^■m2 (g?^i l^g^' (Si) [sm^ ^ (^sm2 [f^';" (S2) [sm'^ 

mi,m'^ .m2 .m2 — — s 

E (2, + i) {^ni,\f^l{Si)\sm[)[s-m,\f^,^{S2)\s-m[) 

nil .m'^ — — s 
^ ^ ( -i\mi— m'j 

= E ~7:^ — ~~r:~ c{sks; m'^qimi) c{sks; -m[qi - ■mi)y/{2ki + l){2k2 + 1) 

^ (2s + 1) 

mi ,m^— — s 

= i-l)>'^+i^6kik2S,i-,2 (4.6) 

where we have made use of Eqs. (2.2), (2.7) and (2.8)). Thus, the P and Q functions assume the simple form 

2s k 

An 



, 2s k 

i'(^i>i;e2,<^2) = ^E(-i)'^'fe E ^^*9(^i' '^1)^^9(^2,02) (4.7) 



fc=0 q=-k 

and 

2s 

47r 



_. 2s k 

O(^i,'/'i;e2,02)-^5](-i)'QL J2Y:q{e,,cj,,)Yk,{e2,h) (4.8) 

fc=0 ij=-fe 

for EPRB spin correlations. 

The F representation can be generalised for pairs of vector statistical variates Xi , X2 denoting 'classical spin 
vectors' constrained by the conditions \Xi\ = ■\/si(si + 1), \X2\ = ■\/s2(s2 + 1) and we obtain, on using the method 
outlined in Section III, 

F{X„X2)= S{^s,{s, + 1) - l^il) E iu I If {^1(^1 + 1)}^^ 

2S2 -, 

X S{^S2{S2 + 1) \X2\) Y. {S2(S2 + 1)}^^ 



k2=^ fe! ■^k2S2 

ki ki 

X E E tql'^KqA9i,MY,:q2i(^2,h) (4.9) 

i?i=-fci gi=-fei 



and the normalised distribution of angular variables i.e., F{9i, <^i; ^2, (^2) is obtained through. 



i^(0i,'^i;02,</'2)= / / xfdXiXldX2F{Xi,x 



Jo 



2si 2s2 ki ki 

= E •^''ifci E •^''^'^^ E E 41'^fe%i(^i>i)^'c;,.(e2>2). (4.10) 

fel— k2—0 qi — — k-i qi—^ki 



For EPRB spin correlations F{di, (j)i;02, (t>2) reduces to 

2s 

47r 



2s k 

F(01,0i;^2>2) = ^E(-l)'-^'fe Y.^kg{0l,MYkq{e2,CJ)2). (4.11) 

Normalisation property of these functions follows naturally from the orthogonality of the spherical harmonics. It 
could be verified that the marginal distributions for system 1 in each of these representations is given by 

P(0i,(/)i) - /dOiP(0i,(/)i;02,'^2) - ^, 

g(0i,(/)i) = /dOig(0i,(/)i;02,02) = ^, 

47r 
i^(^i, (/)i) = Jdn.Fie,, (/-i; 02, 02) - ^, (4.12) 

where df^i = sin^i d9i d(j)i. Similarly, for system 2, the marginal distributions are realised to be P{02,4>2) = 

- — ; Q{02,4'2) ~ -, — ; and F{92,(t>2) = -. — . These marginal distributions are spherically symmetric and hence corre- 

Att Ait An 

spond to totally random orientations (^i, 0i) and (62, 4)2) of the spin vectors associated with the system. 



The spin correlations ( (Si •a)(S2 -b)) could be evaluated using the above distribution functions as follows: We 
make use of the relations between irreducible tensors of rank 1 and the spin components Si, S2, S2, given by 



^i=\/^^^ (fl,(l)-fi(S) 



S2^^^'-^^ (rii(t)+ri(g)) 



^3= f-^^f^$) (4.13) 

which lead to the correspondence rules (see Eqs. (3.34-3.36) of Section III) for the spin operator S in the three 
representations discussed above, as 

P — representation : S — > s n{'K — 9,(j)), 

Q — representation : S — > {s + 1) n(7r — 0, (p), 

F — representation : S — > ys^s + l) n{d, (f)) (4.14) 

where n denotes a unit vector specified by the polar co-ordinates in the arguments. Thus, we have 

((§1 • a)$2 ■ b)) = s^JJdni dn2 P{0i,cl>i;e2,<l>2) (ni • 5) (^2 • &) 

= (5 + 1)2 JJdnidn2 Q{Oi,^i;02, (t>2) (ni • a) (n2 • &) 

= {s(s + 1)}2 // dnidn2 F{0i,4>r,02, h) (m • 5) (^2 • &) (4.i5) 

1 
which, on expressing n ■ a ~ ^ N, Yiq{^)Yiq{'^) ^^^ oii using the orthonormality property of spherical harmonics, 

yeilds the well-known result for EPRB spin-s correlations, 

(Si • a)(S2 • b)) = -^il±ll cos0ab. (4.16) 



Here, cos 0ab — a ■ b. 

One can make use of the addition theorem for spherical harmonics [11] 



Y,yk,{0i,<l>i)Y:^{02,cl>2) - ^^ Pk{cose,2), (4.17) 

Q 

where P^ denote Legendre polynomial of order k; cos0i2 = cos^i cos 6*2 + sin^i sin 02 cos(0i — 02), so that the spin 
distributions reduce to 

1 ''' 
Pi9i2) = ^^E(-l)' (2fc + l) V',, Pfe(cOS0i2) 

1 ^' 
Qi0i2) = (4^E(-1)' (2fc+l) Q?fe Pk{cos9,,) 

1 ^^ 

F(0i2) = 74^2 E("l)' (2fc + l) ^ffe Pfc(cos0i2). (4.18) 

We give below the P, Q, and P functions explicitly for the spin-i case: 

P^(^12)-7^(1-9COS012) 
QH0i2)^ 7^(1-COS012) 

/*(^i2)-^(l-3cos0i2) (4.19) 

where cos0i2 — cos^i cos02 + sin^i sin 02 cos((/)i — ^2). In Fig.l we have plotted P, Q, and F functions for spin 
values s = i, 1, |, 2. It could be observed that the distributions show prominent peaks around 612 — 180° indicating 
the anticorrelation property of associated classical vectors ni and ^2 • 

2s-l 

Expressing (2s ± fc)! = {2s)'^^^'^ TT ( 1 ^ :^ ) etc., in Vsk, Qsfcand Tgk it can be readily observed that 

n—^k 

lim -P.fc = 1, lim Qsk^l, Hm J^.fc = 1, (4.20) 

which, together with the completeness property[ll] 

oc k 

E E YkgidA) Ykq{e\^') = 6{cose-cose') 6(4, ~c^'), (4.21) 

k=Oq=-k 

and the symmetry 

Yk,{e, 0) = {-l)''Yk,{7: -e,7: + (b) (4.22) 

of the spherical harmonic functions, leads to the P, Q, and F distribution functions in the classical limit as, 

lim P(0i,0i;6'2,(/'2) = ^<5(cos6'i+cos6l2) (5(0i - (02 + tt)), 

s—>-oo ^'^ 

lim Q(0i,0i;6'2,02) = ^.^(cos^i + C0S6I2) (5(0i - (02 + tt)), 

s—>-oo ^'^ 

lim F(0i,0i;02,02)= 4^(5(cos0i + cos02)'5(0i-(02+7r)). (4.23) 

Observe that in the classical limit, these distribution functions reflect the perfect anticorrelation between the classical 
vectors ni and n2- 
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FIG. 1. Spin distribution functions (a) P{0\2), (b) Q(6i2) and (c) F(6i2) as a function of the angle 6^12 between the classical 
spin vectors ni and n2 constituting the spin singlet. 
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